Our goals are to compare the CS, N2H + and dust distri- 
butions in a representative sample of high-mass star forming 
dense cores and to determine the physical and chemical prop- 
erties of these cores. We compare the results of CS(5^4) and 
1.2 mm continuum mapping of twelve dense cores from the 
southern hemisphere presented in this work, in combination 
with our previous N2H + (l-0) and CS(2-1) data. We use nu- 
merical modeling of molecular excitation to estimate physical 
parameters of the cores. Most of the maps have several emis- 
sion peaks (clumps). Mean sizes of 17 clumps, having coun- 
terparts in continuum and CS, are 0.30(0.06) pc (continuum) 
and 0.51(0.07) pc (CS). For the clumps with IRAS sources 
we derived dust temperatures: 24-35 K, masses: 90-6900 M Q , 
\Q molecular hydrogen column densities: (0.7 - 12.0) x 10 23 cirT 2 
O and luminosities: (0.6 - 46.0) X 10 4 L Q . LVG densities to- 
, wards CS peaks within the 50" beam (0.56 pc at 2.3 kpc, the 
' average distance of our sample source) vary from source to 
source in the range: (3^40)x 10 5 cirT 3 . Masses calculated from 
' LVG densities are higher than CS virial masses and masses 
■^-s , derived from continuum data, implying small-scale clumpi- 
_ ' ness of the cores. The molecular abundances towards IRAS 
sources in eight objects are X(CS)=(0.3 - 2.7) x 10~ 9 and 
I ■ X(N 2 H + )=(0.3-4.4)xlO- 10 .TheCS and continuum maps have 
been compared with each other and with the N2H + (l-0) maps. 
j ■ For most of the objects, the CS and continuum peaks are close 
r*=* 5 to the IRAS point source positions. The CS(5-4) intensities 
OO correlate with continuum fluxes per beam in all cases, but only 
, in five cases with the N2H + (l-0) intensities. The study of spa- 



■ tial variations of molecular integrated intensity ratios to contin- 
uum fluxes per beam reveals that 7(N2H + )/Fi 2 ratios drop to- 
*G • wards the CS peaks for most of the sources, which can be due to 
£-1 a N 2 H + abundance decrease. For CS(5-4), the I(CS)/F h2 ratios 
Q show no clear trends with distance from the CS peaks, while for 
,H \ CS(2-1) such ratios drop towards these peaks. Possible expla- 
& nations of these results are considered. The analysis of normal- 
. . ized velocity differences between CS and N2FF lines has not 
revealed indications of systematic motions towards CS peaks. 



X 
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Abstract. 



1. Introduction 

The first step in most studies of dense cores of molecular 
clouds is their identification, either by a visual inspection or by 
rather sophisticated numerical analysis of molecular line maps. 
Yet, the maps of "traditional" tracers of dense gas (CS, HCN, 
HCO + , NH3, N2H + etc.) appear to be different from each other 
in many cases. Most distinct differences are often observed be- 
tween CS and N2H + maps. In low-mass starless cores (e.g. 
Tafalla et al. 2002) the CS emission vanishes towards the core 
center while the N2H + emission is still high in the interior re- 
gions. Some more turbulent starless clumps show strong CS 
emission in comparison with N2H + (e.g. Williams & Myers 
1999, Olmi et al. 2005). The optically thick CS lines in the lat- 
ter case are usually asymmetric and blueshifted, implying that 
they are probably tracing infall motions. Both these cases can 
be explained by time dependent chemical models including de- 
pletion molecular species onto grains (Bergin & Langer 1997, 
Bergin et al. 1997, Li et al. 2002, Shematovich et al. 2003, 
Aikawa et al. 2003), suggesting different evolutionary stages 
of clumps within parent molecular cloud. In the cores with em- 
bedded stellar objects, especially those of high mass, the mor- 
phologies of molecular maps become more complicated. This 
is due to the overlapping of the emission coming from the qui- 
escent gas with that from outflows and hot cores, where dust 
grain mantles evaporate because of their proximity to YSOs. 
The N2H + intensities towards these objects may drop, while 
CS peaks towards YSO positions (e.g. Ungerechts et al. 1997, 
Bottinelli & Williams 2004). Thus, chemical models of star- 
less clumps cannot be used to explain the observed differences 
between CS and N2H + maps in high-mass star forming regions. 
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In the last decade, a large sample of high-mass star form- 
ing cores associated with water masers have been mapped in 
the CS(2-1) line (Zinchenko et al. 1994, 1995, 1998). Many of 
these cores have subsequently been mapped in different molec- 
ular lines, including N2H + (l-0) (Pirogov et al. 2003, here- 
after Paper I). Thus, the maps can be compared and important 
quantities can be determined for these regions, such as den- 
sity distribution and chemical composition. Yet, in contrast to 
the N2H + (l-0) lines, which are optically thin in most cases 
(Paper I), the CS(2-1) lines are likely optically thick and possi- 
ble differences between maps in these lines could be connected 
both with optical depth effects and chemical differences. 

In this paper we present the results of 1.2 mm dust con- 
tinuum and CS(5^1) observations in several southern sources. 
Optically thin dust emission is known to trace closely the 
gas component in high density regions and it is insensitive to 
molecular abundance variations and radiative transfer effects. 
This can be used to separate optical depth and chemical effects 
in molecular line maps. 

The goal of the paper is to analyze the dust continuum and 
CS(5-4) data in comparison with our previous CS(2-1) and 
N2H + (l-0) results in order to get reliable information on the 
density and chemical structure of the cores. The paper also con- 
tains estimates of physical parameters including sizes, masses 
as well as density distributions derived from CS(5^1) and 
CS(2-1) data and molecular abundances for several sources. 

2. Observations 

2.1. Source list 

The sources were originally selected from the sample of high- 
mass star forming regions associated with water masers lo- 
cated in the southern hemisphere (Zinchenko et al. 1995, Juvela 
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1996) and from the dense core database of Jijina et al. (1999), 
according to the following criteria: presence of embedded clus- 
ters of stars detected in infrared, and distances to the objects not 
greater than 5 kpc. If no data on clusters were available, we se- 
lected sources with high IR-luminosities (L > 10 4 L ), which 
can be an indirect indication of cluster (Jijina et al. 1999). In 
total, 14 sources have been selected for CS(5^1) observations, 
twelve of them have been observed in continuum at 1.2 mm 
wavelength. Most of the sources have associated IRAS point 
sources. The source list with coordinates and distances is given 
in Tabled 



Table 1. Source list 



Source RA (2000) Dec (2000) D 

(h) r) (s) ( o , „ } (kpc) 



G 264.28+1.48 


08 


56 


27.8 


-43 05 05 


1.4" 


G 265.14+1.45 


08 


59 


24.7 


-43 45 22 


1.7" 


G 267.94-1.06 


08 


59 


03.6 


-47 30 47 


0.7° 


G 268.42-0.85 


09 


01 


54.3 


-47 43 59 


1.3° 


G 269.11-1.12 


09 


03 


32.8 


-48 28 39 


2.6° 


G 270.26+0.83 


09 


16 


43.3 


-47 56 36 


2.6° 


G 285.26-0.05 


10 


31 


30.0 


-58 02 07 


4.7" 


G 291.27-0.71 


11 


11 


49.9 


-61 18 14 


2.7* 


G 294.97-1.73 


11 


39 


12.6 


-63 28 47 


1.2" 


G 305.36+0.15 


13 


12 


33.9 


-62 37 38 


4.2° 


G 316.77-0.02 


14 


44 


58.9 


-59 48 29 


3.1 c 


G 345.01 + 1.80 


16 


56 


45.3 


-40 14 03 


2.1 c 


G 345.41-0.94 


17 


09 


33.7 


-41 35 52 


2.8 r 


G 351.41+0.64 


17 


20 


53.4 


-35 47 00 


L7 d 



Zinchenko et al. (19951 . * Brand & Blitz {1993}, c Juvela (19961 . 
d Neckel (T9781 



2.2. The CS(5-4) observations 

Observations of the CS(5^1) line at 244.9 GHz towards 14 
southern sources have been carried out with the 15-m SEST 
antenna on February, 2001. The N 2 H + (l-0) line at 93.2 GHz 
was observed simultaneously and the results on this line are re- 
ported in Paper I. The telescope half power beam width at the 
CS(5^1) frequency is about 22", and the main beam efficiency 
is 0.5. The system temperature was ~ 400 - 1500 K, depend- 
ing on source elevation and weather conditions. Spectral analy- 
sis was done using an acousto-optical spectrum analyzer (2000 
channels) splitted into two halves to measure the N2H + (l-0) 
and the CS(5^1) lines simultaneously. The frequency resolu- 
tion was 42.6 kHz which corresponds to a velocity resolution 
of 0.052 km s -1 at the CS(5^I) frequency. 

Pointing was regularly checked by SiO maser observations 
and typically it was better than 5". Mapping has been done 
with 20" grid spacing. The data processing included baseline 
subtraction (low-order polynomials) and Gaussian fitting. 



2.3. Observations of dust continuum emission 

In June 2003, twelve sample sources have been mapped in 
dust continuum emission at 250 GHz with the SEST antenna. 
G267.94 and G305.36 were excluded from the list of contin- 
uum observations because they were not completed in CS(5^1) 
andN 2 H + (l-0). 

We used the 37 elements SIMBA bolometer array. The 
HPBW of a single element is about 24", the separation between 
elements on the sky is 44". Mapping has been done by scan- 
ning in azimuth, sampling the signal every 8". The scans were 
spaced by 8" in elevation. This gave the pixel size 8". Typical 
map sizes are 900" x 1200". The raw data have been converted 
by the simbaread software and the processing has been done 
using the MOPSI package written by R. Zylka according to 
the instructions of the SIMBA Observer's Handbook (2003) 
The final r.m.s. noise level derived from regions on the maps 
without sources is 30-60 mJy beam -1 for 9 sources. For three 
sources (G345.01, G345.41 and G351.41) the r.m.s. noise level 
is 160 - 180 mJy beam -1 . The flux measurements in the final 
coadded map are believed to be correct within 20%. Note that 
the continuum fluxes also include contribution from molecu- 
lar lines lying within bolometer bandwidth (mainly CO(2-l)). 
Using the Braine et al. (1995) analysis combined with CO(l- 
0) intensities (Zinchenko et al. 1995), we estimate that such a 
contribution is not significant for any source in our sample. 

3. Results 

3.1. Maps and line parameters 

The CS(5^) integrated intensity maps together with contin- 
uum maps are shown in Fig.^ For comparison, the N2H + (l-0) 
maps (Paper I) are also plotted. In general, the continuum and 
the CS(5-4) maps have similar structures and their morpholo- 
gies range from close to spherical-symmetry (e.g. G270.26) to 
complex clumpy structures with several emission peaks (e.g. 
G3 16.77). IRAS sources are indicated by stars and the uncer- 
tainty ellipses corresponding to 95% confidence level in posi- 
tional accuracy are also shown. For most of the objects (except 
G269.1 1, G316.77 and, probably, G265.14 and G345.41) these 
ellipses overlap with the CS and dust 90% intensity contours, 
indicating coincidence of CS, continuum peaks and IRAS point 
source positions. 

Offset coordinates of CS(5^4) integrated intensity peaks, 
the values of integrated intensities towards these positions and 
parameters of Gaussian fits (main beam temperature, velocity 
and line width) are given in Table|2]with the corresponding 1 <x 
errors, as defined by the fits, in brackets. 

3.2. Clump parameters 

The CS(5-4) and continuum maps have been deconvolved into 
individual clumps using our 2D Gaussian fitting program and 
the method described in Paper I. No attempts to separate indi- 
vidual clumps have been made in G264.28 where local emis- 
sion peaks of nearly equal intensities are located too close to 

1 http://puppis.ls.eso.org/stan7simba/manual/simba/index.html 
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10-1 I -1 

Fig. 1. Continuum 1.2 mm maps (in color) together with CS(5^1) (solid blue contours) and N2H + (l-0) (dashed yellow contours, 
taken from Paper I) integrated intensity maps. Intensity contours range from 10% to 90% of peak values with 10% step plus 5% 
contour in addition (continuum), from 20% to 90% of peak values with 10% step (CS), and from 30% to 90% of peak values with 
20% step (N2H + ). The CS peak integrated intensities are given in Table|2] The peak continuum fluxes (in Jy beam -1 ) are: 1.46 
(G264.28), 2.32 (G265.14), 7.43 (G268.42), 3.44 (G269.ll), 3.34 (G270.26), 4.39 (G285.26), 18.83 (G291.27), 1.64 (G294.97), 
3.37 (G3 16.77), 3.93 (G345.01), 3.39 (G345.41), 12.54 (G351.41). The continuum clumps are marked by numbers as in Tableg| 
IRAS point sources are marked by red stars. The uncertainty ellipses corresponding to 95% confidence level in IRAS point source 
position are also shown. The SEST beam size at the CS(5^4) frequency is shown on the G 264.28 map 



each other. We have not processed the CS(5^4) map in G267.94 
as it has not been completed. The parameters of individual 
CS(5-4) clumps, including relative coordinates of clump cen- 
ters, aspect ratios (the ratios of the extents of the fitted ellip- 
tical Gaussians) and deconvolved angular (A0) and linear (d) 
sizes estimated as geometric mean of the extents of the ellip- 



tical Gaussians at half maximum intensity level are given in 
Table [5] (columns 2-6). The clumps which belong to the same 
object are marked by numbers. 

In columns 7 and 8 of Table|3]the CS(5-4) mean line widths 
and virial masses are given. For homogeneous spherically- 
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Fig. 1. continued 



symmetric clumps with no external pressure and no magnetic 
field, virial masses are given by: 



M vir (M ) = 105 (AV) 2 ■ d , 

where (AV) is the CS(5^1) mean line width (in km s _1 ), de- 
fined as the weighted average of line widths at different posi- 
tions within half-maximum intensity region; d is the CS(5^)-) 
emission region size in pc. Virial masses calculated according 
to the above formula lie in the range ~ 140- 1630 M G , with a 
mean value of 753 M Q . We point out that these values can be 
overestimated, given that the CS(5^4) lines can be broadened 
due to optical depth. Shirley et al. (2003) who have studied a 
large sample of high-mass star forming regions in the northern 
hemisphere found that CS(5-4) line widths are on average 1.3 
times higher than those of presumably optically thin C 34 S(5-4) 
lines. In addition, if density in the cores decreases outwards, 
virial mass should be multiplied by a factor ^^-p) ' wnere P i s 
the power-law index for density radial profile. For p close to 2 
(Paper I) this factor is close to 0.6. 



The columns 2-6 of Table 0] contain parameters of contin- 
uum clumps similar to those of CS clumps. Total fluxes (f total) 
are given in column 7. A presence or an absence of an IRAS 
point source within the half maximum intensity level of indi- 
vidual clumps is marked by plus or minus signs in the last col- 
umn of Tableland Table0] 

All the CS(5^1) clumps have corresponding counterparts 
in the continuum. Their centers coincide within 10" in most 
cases. Larger angular distances (up to 16") between CS and 
continuum clump centers have been found in G294.97(2) and 
G316.77(l), however, in these cases the sensitivity of the CS 
data is rather low. No attempts to reveal individual clumps have 
been done in the case of G264.28 (both for continuum and CS 
maps) and G291.27 (forCS map) where the emission peaks are 
too close to each other and have comparable intensities. 

In several cases (G265. 14, G285.26, G291.27, G345.01 and 
G351.41) the continuum maps show more detailes than their 
CS(5^)-) counterparts probably due to a higher signal-to-noise 
ratio, better spatial sampling and absence of optical depth ef- 
fects. In G351.41, the continuum map reveals two clumps sep- 
arated by about 110" in the north-south direction which are 
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Table 2. The CS(5-4) line parameters 



Source (Aa", AS") I T M b Vlsr AV 

(Kkms- 1 ) (K) (kms- 1 ) (km s" 1 ) 



G 264.28+1.48 


(0,-20) 


13.3(0.2) 


5.08(0.06) 


5.93(0.01) 


2.34(0.03) 


G 265.14+1.45 


(-20,0) 


20.9(0.6) 


7.34(0.18) 


7.80(0.03) 


2.58(0.08) 


G 267.94-1.06 


(0,20) 


16.8(1.1) 


4.01(0.21) 


3.33(0.11) 


4.37(0.27) 


G 268.42-0.85 


(0,-20) 


30.8(0.7) 


10.26(0.22) 


3.51(0.03) 


2.90(0.07) 


G 269.11-1.12 


(0,40) 


35.8(0.7) 


5.82(0.11) 


9.74(0.05) 


5.78(0.13) 


G 270.26+0.83 


(-20,20) 


25.3(0.7) 


6.01(0.17) 


9.67(0.06) 


4.17(0.14) 


G 285.26-0.05 


(0,-20) 


31.8(0.3) 


5.76(0.05) 


3.67(0.02) 


5.04(0.05) 


G 291.27-0.71 


(-20,-40) 


43.0(0.6) 


10.44(0.14) 


-23.75(0.03) 


4.08(0.07) 


G 294.97-1.73 


(20,-20) 


14.9(0.5) 


4.71(0.17) 


-8.01(0.05) 


2.78(0.12) 


G 305.36+0.15 


(20,20) 


18.4(0.5) 


3.66(0.08) 


-38.56(0.05) 


5.46(0.14) 


G 316.77-0.02 


(220,-60) 


22.0(0.3) 


6.10(0.08) 


-38.44(0.02) 


3.41(0.05) 


G 345.01 + 1.80 


(20,-20) 


44.1(0.5) 


8.03(0.08) 


-13.48(0.03) 


5.32(0.06) 


G 345.41-0.94 


(20,0) 


69.2(0.5) 


13.87(0.10) 


-21.16(0.02) 


4.82(0.04) 


G 351.41+0.64 


(0,0) 


116.8(0.6) 


17.19(0.08) 


-7.03(0.01) 


6.42(0.04) 



Table 3. Physical parameters of the CS(5^1) clumps 



Source Aa AS Aspect A0 d (AV) M vil IRAS 

(") (") ratio (") (pc) (kms- 1 ) (M ) 



G 264.28+1.48 
G 265.14+1.45 (1) 
G 265.14+1.45 (2) 
G 268.42-0.85 
G 269.11-1.12 
G 270.26+0.83 
G 285.26-0.05 (1) 
G 285.26-0.05 (2) 
G 291.27-0.71 
G 294.97-1.73 (1) 
G 294.97-1.73 (2) 
G 305.36+0.15 
G 316.77-0.02(1) 
G 316.77-0.02 (2) 
G 316.77-0.02 (3) 
G 316.77-0.02 (4) 
G 316.77-0.02 (5) 
G 345.01 + 1.80 
G 345.41-0.94 
G 351.41+0.64 



4(3) 


-26(1) 


-29(1) 


0(1) 


12(2) 


-96(1) 


-6(1) 


-12(2) 


4(2) 


35(2) 


-15(2) 


20(1) 


-7(1) 


-20(1) 


-49(2) 


-58(1) 


-8(2) 


-3KD 


16(3) 


-13(3) 


54(7) 


0(3) 


15(2) 


33(3) 


216(1) 


-59(1) 


156(2) 


-64(2) 


93(4) 


-17(5) 


17(3) 


-12(2) 


-14(3) 


38(5) 


22(2) 


-23(1) 


4(2) 


-2(1) 


-3(1) 


-5(1) 



3.2(0.8) 42(4) 

1.8(0.1) 58(2) 

1.6(0.2) 57(4) 

1.6(0.1) 64(3) 

1.2(0.2) 46(4) 

1.5(0.2) 38(3) 

1.2(0.2) 18(2) 

1.9(0.3) 34(3) 

1.4(0.1) 70(2) 

1.3(0.3) 51(5) 

1.6(0.5) 63(9) 

1.6(0.5) 37(5) 

1.1(0.2) 32(3) 

2.7(0.6) 32(4) 

2.1(0.6) 62(9) 

1.5(0.3) 57(6) 

2.2(1.0) 39(8) 

1.9(0.1) 85(2) 

1.8(0.2) 59(3) 

1.1(0.1) 35(2) 



0.28(0.03) 2.2(0. 

0.48(0.02) 2.2(0. 

0.47(0.03) 1.7(0. 

0.40(0.02) 2.7(0. 

0.58(0.05) 2.9(0. 

0.48(0.04) 2.6(0. 

0.42(0.04) 4.5(0. 

0.78(0.06) 3.4(0. 

0.92(0.03) 4.0(0. 

0.29(0.03) 2.2(0. 

0.37(0.06) 2.4(0. 

0.76(0.11) 3.7(0. 

0.48(0.04) 3.9(0. 

0.49(0.06) 4.4(0. 

0.93(0.14) 4.1(0. 

0.85(0.09) 3.3(0. 

0.59(0.13) 3.7(0. 

0.87(0.02) 4.0(0. 

0.81(0.04) 3.6(0. 

0.29(0.02) 5.2(0. 



2) 


142 


+ 


2) 


238 


+ 


1) 


141 




2) 


307 


+ 


3) 


506 




3) 


337 


+ 


3) 


868 


+ 


2) 


953 




1) 


1550 


+ 


2) 


152 


+ 


2) 


219 




3) 


1090 




3) 


762 


+ 


2) 


973 


+ 


4) 


1630 




3) 


985 




3) 


870 




1) 


1450 


+ 


2) 


1070 


+ 


1) 


814 


+ 



associated with the well-known massive star forming regions 
NGC 6334 I and NGC 6334 I(N), thought to be at different 
evolutionary stages (e.g. McCutcheon et al. 2000). Although 
the northern clump in our CS(5-4) map is not so prominent 
as in continuum, the map shows similar morphology and is in 
general agreement with the CS(7-6) results from McCutcheon 
et al. (2000). 

Angular and linear sizes of 20 CS clumps lie in the ranges: 
A©= 18-85" andc/ = 0.28 -0.93 pc with <A0) = 50"(6") and 
(d) = 0.52(0.07) pc; the numbers in brackets are r.m.s. devia- 
tions from the mean. The sizes of 25 nearby continuum clumps 
are: A0 = 12 -56" and d = 0.13 -0.84 pc with (A0) = 33(4)" 



and (d) = 0.34(0.06) pc. The continuum clumps in most of 
the sources are smaller than the CS clumps. For 17 clumps de- 
tected both in the continuum and CS (excluding G264.28 and 
G291.27) {d) = 0.30(0.06) pc (continuum) and 0.51(0.07) pc 
(CS). The difference in sizes could be connected with low 
signal-to-noise ratios and undersampling of the CS maps which 
may produce enlarged CS clumps sizes. Higher signal-to-noise 
ratios in the continuum observations allowed us to reveal mul- 
tiple clumps in some sources (e.g. G269.1 1, G291.27) while in 
CS we see single large clump. Yet, it is also possible that the 
differences in sizes could be due to the CS opacity or abun- 
dance gradients within clumps. The CS clump sizes are close 
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Source 


Aa 


AS 


Aspect 


A0 


d 


F total 


IRAS 




(") 


(") 


ratio 


(") 


(PC) 


(Jy) 




G 264.28+1.48 


2.2(0.4) 


-33.0(0.2) 


4.0(0.2) 


36(1) 


0.24(0.01) 


5.0 


+ 


G 265.14+1.45 (1) 


-28.1(0.3) 


3.2(0.3) 


1.4(0.1) 


30(1) 


0.25(0.01) 


5.1 


+ 


G 265.14+1.45 (2) 


14.4(0.6) 


-88.7(0.6) 


1.3(0.1) 


50(1) 


0.42(0.01) 


5.8 




G 265.14+1.45 (3) 


-8.1(0.3) 


-38.6(0.5) 


1.5(0.1) 


36(1) 


0.30(0.01) 


5.0 


+ 


G 268.42-0.85 


-5.9(0.2) 


-6.9(0.2) 


1.37(0.03) 


38.2(0.4) 


0.241(0.003) 


24.8 


+ 


G269.ll-l.12 (1) 


1.1(0.2) 


34.5(0.2) 


1.6(0.1) 


27(1) 


0.34(0.01) 


7.5 




G 269.11-1.12 (2) 


24(1) 


82(1) 


1.4(0.1) 


44(2) 


0.56(0.03) 


2.9 




G 270.26+0.83 


-21.6(0.2) 


20.9(0.2) 


1.2(0.1) 


20.8(0.4) 


0.26(0.01) 


6.6 


+ 


G 285.26-0.05 (1) 


-6.7(0.1) 


-11.3(0.1) 


1.03(0.03) 


19.8(0.3) 


0.45(0.01) 


7.0 


+ 


G 285.26-0.05 (2) 


^18.3(0.4) 


-56.7(0.3) 


1.5(0.4) 


12(2) 


0.27(0.03) 


1.7 




G 291.27-0.71 (1) 


11.0(0.1) 


-36.2(0.1) 


1.28(0.02) 


41.9(0.3) 


0.548(0.004) 


74.5 


+ 


G 291.27-0.71 (2) 


-74(1) 


-92(1) 


2.5(0.2) 


45(2) 


0.59(0.03) 


6.2 




G 291.27-0.71 (3) 


-121(1) 


-183(1) 


1.1(0.1) 


50(1) 


0.65(0.02) 


5.0 




G 294.97-1.73 (1) 


6.4(0.4) 


-12.5(0.4) 


1.4(0.1) 


44(1) 


0.254(0.004) 


5.9 


+ 


G 294.97-1.73 (2) 


65(1) 


15.7(0.4) 


1.1(0.1) 


42(1) 


0.24(0.01) 


3.9 




G 316.77-0.02 (1) 


200(1) 


-61.8(0.3) 


2.7(0.1) 


46(1) 


0.69(0.01) 


9.1 


+ 


G 316.77-0.02 (2) 


168(1) 


-66.2(0.4) 


2.5(0.2) 


41(1) 


0.62(0.02) 


8.7 


+ 


G 316.77-0.02 (3) 


95(1) 


-24(1) 


1.8(0.1) 


56(1) 


0.84(0.02) 


6.2 




G 316.77-0.02 (4) 


21(1) 


-15(1) 


1.1(0.1) 


40(1) 


0.60(0.02) 


5.7 


- 


G 316.77-0.02 (5) 


-14(1) 


25(1) 


1.3(0.1) 


39(2) 


0.59(0.02) 


3.6 




G345.01 + 1.80 (l) 


15(1) 


-22.0(0.4) 


2.4(0.1) 


42(1) 


0.43(0.01) 


14.3 


+ 


G345.01 + 1.80(2) 


-69(1) 


22(1) 


1.9(0.2) 


48(2) 


0.49(0.02) 


6.6 




G 345.41-0.94 


10(1) 


-7.7(0.4) 


2.3(0.1) 


42(1) 


0.57(0.02) 


12.6 


+ 


G 351.41+0.64 (1) 


-5.2(0.3) 


-5.5(0.2) 


2.1(0.2) 


16(1) 


0.13(0.01) 


19.9 


+ 


G351.41+0.64 (2) 


20.3(0.3) 


99.8(0.3) 


1.05(0.04) 


32(1) 


0.26(0.01) 


21.8 





to those of N2H + ones (Paper I). Aspect ratios for CS and con- 
tinuum clumps are: 1.1 - 3.2 and 1.0 - 4.0, with mean val- 
ues 1.5(0.1) and 1.3(0.2), respectively. This parameter could 
be used for selecting nearly circular cores and exploring radial 
dependences of physical parameters within the clumps. Such 
an analysis lies beyond the scope of the present paper and we 
postpone it to future publications. 

4. Physical parameters derived from continuum 
data and from LVG modeling 

Using our 1.2 mm data and IRAS fluxes it is possible to es- 
timate physical parameters of the clumps including dust tem- 
peratures, masses, hydrogen column densities and luminosities. 
The CS(5-4) data together with the (2-1) data (Zinchenko et al. 
1995) allowed us to calculate LVG densities, CS column den- 
sities and molecular abundances for eight sources. 

4.1. Dust temperatures, masses, hydrogen column 
densities and luminosities 

For the clumps with nearby IRAS sources we have derived 
dust color temperatures by fitting the IRAS and 1.2 mm total 
fluxes (f tota i, Table0} with a 2-temperature greybody curve and 
assuming optically thin conditions: A[V^B v (Ti) + A2V^B V (T2) 
(Mozurkewich et al. 1986). We exclude G269. 1 1 from the anal- 
ysis because the IRAS source lies outside the half-maximum 
level in this source. We set the power-law index of dust 



emissivity-frequency dependence (B) equal to 2 in our calcu- 
lations. Setting B — 1 leads to 10-20 % increase in dust temper- 
atures. The values of dust temperatures (Tj, cold component) 
are given in column 2, Table [5] The temperatures of the hot 
dust component lie in the range: ~ 95 - 115 K. Dust masses 
have been calculated according to the expression (e.g. Doty & 
Leung 1994): 

. . ^total D 2 

M d = , 

kx.2Bx.2iTi) 

where D, k\ i and fli.2(7d) are the source distance, the dust mass 
absorption coefficient and Plank function at 1.2 mm, respec- 
tively. Ossenkopf & Henning (1994) derived a mass absorption 
coefficient of 1 cm 2 g -1 at 1.3 mm wavelength appropriate for 
cold dust grains covered with thick icy mantles. We adopted 
this value for our calculations assuming that water ice cannot 
evaporate at the derived values of Td (e.g. van Dishoeck 2004). 
Gas masses (M g ), calculated assuming a gas-to-dust mass ratio, 
R,„ = 100, lie in the range -90-6900 M Q . Note thatthese val- 
ues could be underestimated if the dust temperature decreases 
outwards as expected in the case of internal heating. Using the 
simple model of optically thin spherical dust shell with radial 
temperature and density gradients (Pirogov & Zinchenko 1998; 
see the Appendix) and taking the fluxes at 1.2 mm and 100 fim 
for each source it is possible to calculate the ratio of masses 
obtained in this model and in the isothermal approximation. It 
mainly depends on the ratio of outer-to-inner radii of the shell 
and the density-radius power-law index, increasing when these 
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two parameters increase. Taking the power-law index p — 2 and 
a radii ratio of 10 4 we find that the mass of the shell for some 
sources can be up to 5 times higher than the mass calculated in 
the isothermal approximation. 

Hydrogen column densities have been calculated from peak 
fluxes per beam and dust temperatures (e.g. Motte et al. 1998): 



nmR-'k^B^iTi)' 

where £2 is beam solid angle, m = 2.33 amu is mean molecular 
mass. The A^h, values given in column 5 of Table[5]fall in the 
range: (0.7- L2.0)x 10 23 cm 4 . 



Table 5. Physical parameters derived from continuum data 



Source 


r d (K) 


Af g (Af ) 


L(L Q ) 


Af H2 (cm 2 ) 


G 264.28+1.48 


32 


93 


6.6 x 10 3 


7.0 x 10 22 


G 265.14+1.45 


30 


624 


5.8 x 10 4 


1.2 x 10 23 


G 268.42-0.85 


35 


360 


4.6 x 10 4 


3.2 x 10 23 


G 270.26+0.83 


29 


469 


1.8 x 10 4 


1.8 x 10 23 


G 285.26-0.05 (1) 


33 


1418 


4.6 x 10 5 


2.0 x 10 23 


G 291.27-0.71 (1) 


25 


6910 


2.9 x 10 5 


1.2 x 10 24 


G 294.97-1.73 (1) 


27 


98 


5.9 x 10 3 


9.6 x 10 22 


G 316.77-0.02 


28 


2415 


1.9 x 10 s 


1.8 x 10 23 


G 345.01 + 1.80 (1) 


30 


639 


7.4 x 10 4 


2.0 x 10 23 


G 345.41-0.94 


29 


1035 


3.4 x 10 5 


1.8 x 10 23 


G 351.41+0.64 (1) 


30 


591 


8.5 x 10 4 


6.5 x 10 23 



The comparison of M s masses with CS (Table|3J and N2H + 
virial masses (Paper I), for closely located clumps, shows that 
they agree within a factor of 2 except for G291.27, where M„ 
exceeds CS virial mass by a factor of 4.5. As the uncertainties 
in dust masses (mainly due to uncertainties in distance and in 
dust absorption coefficient and due to deviations from isother- 
mal approximation) can be rather high, we conclude that the 
two mass estimates are in reasonable agreement. 

Bolometric luminosities (L) calculated by integrating the 
fitting curves over total frequency range are given in column 4 
of Table0 They lie in the range: (0.6 - 46.0)x 10 4 L indicating 
the presence of embedded high-mass objects. The L/M g ratio 
which is considered to be proportional to star formation rate 
ranges from 38 to 325 L G /M Q for the sources in Table [5] with 
a mean value 139 + 105 L Q /M G . This is consistent with (L/M) 
estimates for large samples of high-mass star forming regions 
with masses obtained from dust continuum data: 71 +56 L e /M Q 
(Faundez et al. 2004), 136 + 100 L /M (Shirley et al. 2003), 
120 + 90 Lq/Mq (the data from Beuther et al. 2002 rescaled by 
Mueller et al. 2002). 

4.2. LVG densities and molecular abundances 

Taking into account peak line intensities of two different CS 
transitions one can derive hydrogen densities for given kinetic 
temperature using the LVG approach. 

Nine sources of our sample were previously mapped in the 
CS(2-1) line with the SEST antenna (Zinchenko et al. 1995) 



with a 50" beam and 40" grid spacing. The peak CS(2-1) po- 
sitions were also observed in the CO(1-0) line by these authors. 
We have compared their data with our CS(5^1) data for eight 
sources (except G267.94). 

In order to use both the CS(2-1) and the CS(5-4) data in 
model calculations it is important to convolve them to the same 
beam. For this purpose we have smoothed the CS(5^4) main 
beam temperature maps to the 50" resolution using our 2D 
Gaussian fitting program and have calculated errors of the con- 
volved CS(5^4) temperatures using the propagation of errors 
approach. We applied the CS-H2 collisional rates from Turner 
et al. (1992) which are tabulated for given kinetic temperatures. 
The model is isothermal and kinetic temperatures have been es- 
timated according to the temperatures of the CO(1-0) line and 
of the dust (cold component) (Table[5}- The assumed values of 
kinetic temperatures are given in Table|6] Twelve CS rotational 
levels have been considered in the calculations. 

The results of model calculations for eight analyzed sources 
are shown in Fig. 13 as density versus projected distance from 
the positions of peak CS(5^4) integrated intensities (Tables- 
Central densities vary in the range: (3 - 40) x 10 5 cirT 3 . In 
all but two sources (G270.26 and G294.97) density clearly de- 
creases with distance from the CS(5-4) peaks. In four sources 
(G265.14, G268.42, G269.ll and G291.27) density falls by 
about an order of magnitude at ~ 0.8-2 pc from the center. The 
slopes of density-radius dependences for the outer regions of 
these sources vary from -2.0 to -3.6. In G265.14 and G291.27 
these trends may be due to the existence of neighboring low 
density clumps (see Fig.[Q. Note that the calculated LVG den- 
sities depend on the adopted kinetic temperatures (the lower 
the temperature, the higher the density). If the kinetic tempera- 
ture decreases outwards, calculated densities far from the cen- 
ter could be underestimated. Changing the temperature from 
40 K to 20 K leads to a density increase by about half an order 
of magnitude. Also, densities could be slightly overestimated 
towards the positions of IRAS sources due to the influence of 
infrared pumping on CS excitation (see Section[6). 

Taking into account density estimates, it is possible to cal- 
culate gas masses for the most intensive CS(5^4) clumps in 
each source. Assuming spherically-symmetric clumps with ra- 
dius d/2, where d is the CS(5^4) clump size (Table [3j, the 
masses could be calculated as: M = 4^/3 m (n) (d/2) 3 , where 
(n) is the density averaged over positions within the area of ra- 
dius d/2. Using these masses we calculate mass ratios, M v i r /M, 
that should be considered as dense gas volume filling factor. 
The mean mass ratio for eight sources is 0.2 + 0.2. This value 
is consistent with the one found by Shirley et al. (2003) for 
their sample of about 40 dense high-mass star forming regions 
(0.5 + 0.7). It is also consistent with the estimate ^ 0.2 for 
dense gas volume filling factor found by Juvela (1998) from 
multi-line CS and C 34 S profile modeling towards a sample of 
southern high-mass star forming regions. 

Using CS column densities from LVG calculations and H2 
column densities from 1.2 mm continuum data smoothed to 
the 50" beam, the CS abundances have been calculated for 
eight sources towards positions close to IRAS sources. The 
calculated values are given in Table [6] For the same posi- 
tions we have calculated the N2H + LTE column densities (at 
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Fig. 2. LVG densities for eight sources versus projected dis- 
tance from the CS(5^) peaks. Error bars correspond to 1 <x un- 
certainty in the CS line temperatures 



?ex = 10 K) and the N2H + abundances. The abundances lie in 
the following ranges: X(CS)=(0.3-2.7)xlO~ 9 , X(N 2 H + )=(0.3- 
4.4) x 1CT 10 . The lowest X(N 2 H + ) values have been found to- 
wards the most luminous IRAS sources. Mean N 2 H + abun- 
dances for several clumps were calculated in Paper I using 
N 2 H + virial mass estimates. There is a reasonable agreement 
between different X(N 2 H + ) estimates for clumps in G265.14, 
G269.1 1 and G294.97. For G268.42 our estimate of X(N 2 H + ) 
is about 4 times lower than the estimate from Paper I. For 
G285.26 and G291.27 our X(N 2 H + ) values are about an order 
of magnitude lower than the estimates from Paper I, but they 
relate to different clumps, implying abundance variations over 
the sources. The role of abundance variations is discussed in 
Section[6] 



Table 6. Molecular column densities and abundances smoothed 
over 50" Gaussian beam towards positions close to IRAS 



Source, 


7kin 


N(H 2 ) 


7V(CS) 


X(CS) 


N(N 2 H + ) 


Z(N 2 H 


position 


(K) 


(cm- 2 ) 


(crrr 2 ) 


(10-») 


(cnr 2 ) 


(io- 10 






(10 22 ) 


(10 13 ) 




(10 12 ) 




G 264.28 (0,-40) 


20 


3.1 


3.4 


1.1 


4.1 


1.3 


G 265.14 (0,-20) 


20 


4.7 


4.8 


1.0 


16.1 


3.4 


G 268.42 (0,0) 


40 


13.8 


14.0 


1.0 


8.7 


0.6 


G269.ll (0,0) 


20 


3.4 


4.5 


1.3 


14.0 


4.2 


G 270.26 (0,0) 


40 


2.9 


7.9 


2.7 


13.0 


4.4 


G 285.26 (0,0) 


40 


6.0 


4.4 


0.7 


1.6 


0.3 


G 291.27 (20,-20) 


40 


54.5 


17.5 


0.3 


15.6 


0.3 


G 294.97 (0,0) 


40 


4.4 


5.2 


1.2 


1 1.0 


2.5 



Note, that other physical parameters found for the sample 
cores (sizes, aspect ratios, densities, masses, dust temperatures 
and luminosities) lie within the ranges found in the large sur- 
veys of high-mass star forming regions performed during re- 
cent years in CS(5-4) (Plume et al. 1997, Shirley et al. 2003) 
and 1.2 mm continuum (Beuther et al. 2002, Mueller et al. 
2002, Faundez et al. 2004, Fontani et al. 2005). 

5. Comparison of molecular and continuum data 

The maps in Fig.^show in many cases significant differences 
between CS(5-4) and N 2 H + (l-0) intensity distributions as well 
as between 1.2 mm continuum and N 2 H + (l-0). We have per- 
formed detailed point-to-point comparison of molecular inten- 
sities and continuum fluxes for each source and the results are 
given below. In addition, we have compared the CS(5^4) and 
the N 2 H + (l-0) line velocities and line widths. 

5.1. CS and dust maps 

The comparison of the CS(5-4) and 1.2 mm continuum data 
shows that in most cases the CS intensities depend almost lin- 
early on continuum fluxes implying that CS(5^1) is a good 
tracer of total gas column density. Some deviations from linear 
dependences could be connected with differences in excitation 
conditions in clumps within the same source. 

5.2. CSandN 2 H + maps 

In many cases (especially in G285.26, G291.27, G345.01, 
G345.41 and G351.41, see Fig. □ the N 2 H + (l-0) and the 
CS(5^1) intensities have different spatial distributions which, 
in general, could be caused by opacity effects, line excitation 
or differences in chemical composition of individual clumps. 

For comparison with the N 2 H + (l-0) data obtained with the 
55" beam, the CS(5-4) integrated intensities have been aver- 
aged over areas with 25" radius around each position. In five 
objects (G265.14, G268.42, G269.ll, G270.26 and G294.97) 
/(CS) almost linearly depends on 7(N 2 H + ). Other sources show 
a low correlation or no correlation. In some of them (G285.26, 
G345.01 and G351.41) the plots split into separate branches 
which probably belong to different clumps. In Fig.[3]the CS(5- 
4) integrated intensities are shown versus 7(N 2 H + ) for three 
representative sources with high correlation (G269.1 1) and low 
correlation (G285.26 and G291.27) between intensities. 

Comparison of the N 2 H + (l-0) maps with the CS(2-1) 
maps (Zinchenko et al. 1995) shows that in some cases 
(G264.28, G285.26 and G291.27) they are in better agreement 
than with the CS(5-4) maps. 

5.3. Intensity ratios 

In order to study spatial variations of molecular and contin- 
uum emission within the sources we have calculated the CS(5- 
4) to the N 2 H + (l-0) integrated intensity ratios as well as the 
ratios of molecular integrated intensities to continuum fluxes 
(/(CS)/Fi. 2 and /(N 2 H + )/Fi. 2 ) for each source position. The 
continuum fluxes (F1.2) have been averaged over areas with 
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Fig. 3. The CS(5^) versus N2H + (l-0) integrated intensities (both in K km s ') for three representative sources. The CS data is 
smoothed to the 50" resolution 



25" radius around each N2H + position for 7(N2H + )/Fi .2 depen- 
dences. The ratios of molecular integrated intensities to contin- 
uum fluxes give information on abundance variations over the 
sources, if lines are optically thin as in the case of N2H + (l-0) 
(Paper I). We used data with signals higher than 3<x, and re- 
jected map positions with intensities below 5% the peak values 
(3% for G268.42 and G291.27), to avoid spikes of ratios which 
sometimes occur at the edges of the maps. 

The dependences of 7(CS)//(N2H + ) (upper panels), 
7(CS)/fi. 2 (middle panels) and I(N 2 U + )/F h2 (lower panels) 
with respect to projected distance from the CS(5-4) peaks (b) 
are shown in Fig. |4] for three selected sources, where differ- 
ences between the CS(5-4) and N2H + (l-0) intensity distribu- 
tions are among the most prominent. In several sources (espe- 
cially in G285.26, G345.41 and G351.41) the 7(CS)//(N 2 H + ) 
ratios have a tendency to decrease with distance. Six sources 
demonstrate variations of 7(CS)//(N2H + ) across the source by 
more than an order of magnitude (G285.26, G291.27, G345.01, 
G345.41 and G351.41). The highest ratios (~ 13-15) are ob- 
served in G285.26, G291.27 and G345.41, in the vicinity of 
the most luminous IRAS sources of our sample. The lowest 
ratios 0.1) are observed in G3 16.77, G345.01 and G351.41 
towards positions close to the centers of the N2I-F clumps with- 
out CS counterparts and IRAS sources. 

There is a rather large scatter both in 7(CS)/Fi 2 and 
7(N2H + )/Fi.2 ratios across the sources (about an order of mag- 
nitude). The /(N2H + )/Fi 2 ratios clearly drop towards CS peaks 
in most cases while the 7(CS)/Fi.2 ratios demonstrate no such 
trends. No clear trends of /(N2H + )/Fi.2 versus b have been 
found in G268.42 and G294.97. 

For several sources we have calculated ratios of the CS(2- 
1) both to the N2H + (l-0) integrated intensities and to contin- 
uum fluxes averaged over areas with 25" radii around each 
CS(2-1) position. It is found that these ratios show similar 
trends with b as the 7(CS)//(N 2 H + ) ratios for CS(5^1) and 
the 7(N2H + )/Fi .2 ratios, respectively (see plots for G285.26 in 
Fig. For CS(2-1) an increase of the 7(CS)/Fi. 2 ratios with 
b is probably connected with line saturation (see Section[6]and 
AppendixlAl. 

The comparison of intensity ratios at positions close to 
IRAS sources with bolometric luminosities (Table |5} shows 
that the 7(CS)/7(N2H + ) ratios are clearly enhanced towards the 



three most luminous sources (G285.26, G291 .27 and G345.41), 
the 7(N2H + )/Fi .2 ratios drop towards these sources. 

5.4. CS and N 2 H + velocities and line widths 

The velocity comparison of two molecular lines with different 
optical depths gives information about internal kinematics of 
the sources. Mardones et al. (1997) have introduced normal- 
ized velocity difference: 6V = (V,hick - V thin) I thin, where 
y thick and Vthin are optically thick and thin line velocities, re- 
spectively, AV,hi„ - optically thin line width. Their analysis of 
molecular lines in low-mass cores has proved this parameter to 
be an effective tracer of infall (Mardones et al. 1997, Lee et al. 
1999, 2001) and/or core accretion motions (e.g. Tafalla et al. 
2002). Although N 2 H + and CS are likely to trace different re- 
gions, as in low mass cores, it is still interesting to quantify the 
relative motions between different parts of the same cloud and 
see if there are signs of systematic contraction or expansion. 

The N2I-F lines most probably are optically thin in the 
sources (Paper I) while the CS lines could be optically thick. 
In order to search for systematic motions we have calculated 
normalized velocity differences between CS and N2H + (5 Vcs) 
taken from Gaussian fits for each source position. Both the 
CS(5-4) and CS(2-1) data (Zinchenko et al. 1995) have been 
used. Given that the signal-to-noise ratios of the CS(5^4) data 
are rather low, we have only analyzed the spectra with inte- 
grated intensities not lower than 20%-30% of the maximum 
value. Even if our spatial resolution does not permit to resolve 
the regions with systematic motions, one could expect to de- 
tect non-zero normalized velocity differences towards positions 
close to density peaks. 

According to the criterion \6V\ > 0.25 from Mardones et al. 
(1997) no indications of systematic motions have been found 
towards CS(5^4) peaks. However, we have found such indica- 
tions towards some other positions in several sources. In partic- 
ular, in G285.26 the 6Vcs values towards the (0",0") position 
are -0.58(0. 15) and -0.73(0. 19) for CS(5-4) and CS(2-1), re- 
spectively, probably implying infall motions. 

Both CS and N2H + line widths are much higher than 
the thermal widths, and the line profiles in most sources are 
nearly symmetric without prominent self-absorption features. 
The CS(2-1) lines are broader than the N 2 H + (l-0) lines in 
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Fig. 4. The CS(5-4)/N2H + (l-0) integrated intensity ratios and the ratios of molecular intensities to continuum fluxes versus 
projected distance from the CS(5^4) peaks for three sources with most prominent differences between CS and N2H + (in lower 
panels, F 1.2 is the flux at 1.2 mm in Jy beam 1 averaged over areas with 25" radius around each N2H + position). In addition, for 
G285.26 the CS(2-l)/N2H + (l-0) and the CS(2-l)/continuum intensity ratios are shown by red triangles. In this case continuum 
fluxes have been averaged over areas with 25" radii around each CS(2-1) position. Error bars correspond to l<x uncertainties of 
the ratios 



most cases (mean line width ratios vary from source to source 
in the range: 0.9-1.7). The CS(5^4) lines are narrower than 
CS(2-1) and their widths are on average close to the N2H + (1- 
0) widths implying that the CS(5-4) optical depths are not high 
(mean line width ratios vary in the range: 0.7-1.3). In some 
sources (G269.ll, G285.26, G291.27, G345.01 and G345.41) 
the CS(5^) line widths are broader towards CS peaks and 
decrease outwards. Similar trends found in Paper I for the 
N2H + (l-0) line widths were explained as an influence of en- 
hanced dynamical activity in the central regions of the cores 
in the vicinity of IRAS sources. In the case of CS(5^4), these 
trends in addition could be connected with enhanced optical 
depth towards CS peaks. 



The clumps where N2H + drops towards CS peaks 
(Section l5~3l have enhanced CS(5-4) mean line widths 4 km 
s _1 , Table [5} implying either enhanced optical depth or some 
dynamical processes which broaden line profiles towards these 
peaks. There is also a weak correlation between AV(CS) and 
luminosities for positions close to IRAS sources. 



6. Discussion 

Here we examine possible reasons for the observed variations 
of the CS/N 2 H + , CS/dust and N 2 H + /dust ratios. 

Opacity effects. The effect of line saturation due to opacity 
leads to a decrease of the relative line intensities towards CS 
and dust emission peaks (where column density is presumably 
the highest one). We have not seen such effect for the CS(5-4) 
data and this can be considered as an implicit argument against 
high optical depth in these lines. On the other hand a decrease 
of CS(2-1) relative intensities towards density peaks in several 
sources could be connected with saturation. Using the C 34 (2-l) 
data for several sample sources (Zinchenko et al. 1995) we have 
found that the optical depth of the CS(2-1) line is higher than 
unity implying a decrease of the CS(2-1) intensity compared 
with the optically thin case. For G285.26 (-40", -40"), t = 2.5 
and the CS(2-1) intensity decreases 2.7 times. The N2H + (1- 
0) lines most probably are not optically thick (Paper I) within 
the 55" beam and it is hard to consider opacity effects as an 
explanation of the observed variations on the scales resolved 
by the beam. 
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Excitation. Molecular line excitation is determined by col- 
lisions and radiative transfer processes. Differences in critical 
densities of the CS(5^4) and N2H + (l-0) lines by about an order 
of magnitude can produce differences in their excitation and in- 
tensity distributions. The analysis given in Appendix lAlshows 
that density, kinetic and dust temperature gradients can cause 
a decrease of the 7(N2H + )/Fi .2 ratios by a factor not more than 
3 in the density range 10 5 - 10 6 ctrT 3 when Tkin and r<j varies 
from 20 K to 40 K. Yet, they cannot be the reason for the ob- 
served drops of these ratios by an order of magnitude or higher, 
as, for example, in G285.26 where high density contrasts are 
not observed (Fig. [3J and in other sample sources. 

In principle, the observed variations could be connected 
with radiative transfer effects caused by pumping from nearby 
luminous infrared sources which can overpopulate upper rota- 
tional levels of the molecule and underpopulate lower levels. 
This effect in our case could make the CS(5^4) line intensi- 
ties higher and the N2H + (l-0) line intensities lower compared 
with the case of no pumping. However, Carroll & Goldsmith 
(1981) showed that for CS the size of the region affected by 
IR pumping via vibrational transitions is ^ 0.03 pc. We have 
done similar analysis for N2H + using parameters of vibrational 
transitions from Botschwina (1984), Owrutsky et al. (1986) 
and Heninger et al. (2003). It is found that the size of the af- 
fected region considering the bending mode is ;S 0.02 pc and 
is several times lower if one considers pumping via transitions 
in the stretching modes. The details of calculations are given 
in the Appendix [B] A possible radiative pumping via rota- 
tional transitions due to the background radiation field, includ- 
ing IR source, is also found to be unimportant for our sources 
(AppendixlBl. 

Chemistry. As both the opacity and molecular line excita- 
tion effects are found to be insufficient, the observed variations 
of molecular intensity to continuum flux ratios could mainly 
be connected with abundance variations of these molecules. 
Moreover, for different clumps within G285.26 and G291.27 
the N2H + abundances are found to vary by an order of magni- 
tude (see Section l4~2l . High CS-to-N2H + intensity ratios previ- 
ously have been found towards some starless clumps, for exam- 
ple, in S68N (Williams & Myers, 1999) and in Perseus (Olmi 
et al. 2005). They can be explained by N2H + underabundance 
which occurs at early phases of core evolution (Bergin et al. 
1997). However, in our cores the high CS-to-N2H + intensity 
ratios are observed in the vicinity of IRAS sources indicating 
that they are evolved regions. The drop of the N2H + intensity 
towards YSOs previously was observed in Orion (Ungerechts 
et al. 1997) and Cepheus A (Bottinelli & Williams 2004). Hot 
core chemical models (e.g. Nomura & Millar 2004) where ul- 
traviolet radiation from embedded stars leads to evaporation of 
grain mantles and CS enhancement, cannot explain why N2H + 
intensities do not correlate with gas column density in these 
regions. 

Recently Lintott et al. (2005) have proposed an alternative 
model where CS to N2H + abundance ratio becomes enhanced 
due to higher collapse rates compared to free-fall regime (ac- 
celerated collapse). As a result, high densities are achieved 
before CS and the molecules responsible for N2H + removal 
are freezing onto grains. Yet, due to gradients of excitation 



temperatures, systematic collapse motions predicted by this 
model probably should give asymmetric or self-absorbed pro- 
files of optically thick lines (e.g. CS(2-1)) which are not ob- 
served in most of the sources. Normalized velocity differences 
(Section 15. 4> have not indicated such motions towards CS 
peaks in our sources according to the Mardones et al. (1997) 
criterion. Yet, if the sources consist of a large number of small 
unresolved clumps with low volume filling factor as implied 
in Section I4~2l the observed CS lines could have low effective 
optical depth and nearly symmetric profiles. In this case sys- 
tematic collapse motions could be detected as additional line 
broadening towards core centers. The enhanced CS(5^) line 
widths have been detected towards intensity peaks in several 
sources (Section l5.4> where the /(N2H + )/Fi.2 ratios are low. 
If the CS lines actually are not optically thick, these enhance- 
ments could imply some dynamical processes (e.g. collapse 
motions) in gas traced by CS. In order to decide whether or not 
the accelerated collapse model could be applied to these cores, 
one needs to conduct detailed radiative transfer calculations in 
comparison with observational data including other molecular 
tracers (e.g. HCO + , NH 3 , SO). 

Another possible explanation of the N2H + abundance drop 
can be related with its dissociative recombination. Earlier it 
was assumed that dissociative recombination of N2H + leads 
primarily to the formation of molecular nitrogen which can 
form N2H + again via reaction with H3 . However, Geppert et 
al. (2004) found that another channel which leads to forma- 
tion of the NH and H molecules dominates. Perhaps this could 
be the reason why the N2H + abundance drops in vicinities of 
luminous YSOs. It is well known that UV radiation penetrates 
deep into dense molecular clouds, probably due to their clumpy 
structure (e.g. Meixner et al. 1992, Meixner & Tielens 1993). 
Therefore, one can expect an enhanced ionization in the neigh- 
borhood of luminous young stars with a high UV flux, on the 
scales comparable to those which are discussed here. However, 
published results on ionization fraction in massive cores do not 
support this view (Bergin et al. 1999). On the other hand, pre- 
liminary estimates of this fraction from our measurements for 
similar sources in the northern sky do show noticeable varia- 
tions within the sources (Zinchenko et al, in preparation). A 
conclusive decision on the role of this factor in the observed 
chemical differentiation needs further investigation including 
chemical modeling. 



7. Conclusions 

In order to get reliable information on the density and chemi- 
cal structure of the dense cores associated with high-mass star 
forming regions twelve objects from southern hemisphere have 
been mapped in the CS(5^4) line and in dust continuum at 
1.2 mm. We have compared CS(5-4) and continuum data with 
each other as well as with the N2H + (l-0) (Paper I) and the 
CS(2-1) (Zinchenko et al. 1995) data. Besides, physical pa- 
rameters of the cores have been derived from the CS(5-4) and 
the continuum data. 

The results can be summarized as follows: 
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1. Most of the maps have several emission peaks (clumps). 
Mean sizes of 17 clumps, having counterparts in continuum 
and CS, are 0.30(0.06) pc (continuum) and 0.51(0.07) pc 
(CS). The CS virial masses lie in the range ~140- 
1630 M . For the clumps with IRAS sources we derived 
dust temperatures: 24-35 K, masses: 90-6900 M Q , molec- 
ular hydrogen column densities: (0.7- 12.0)x 10 23 ctrT 2 
and luminosities: (0.6-46. 0)x 10 4 L . 

2. Using the CS(5^1) and the CS(2-1) data (Zinchenko et 
al. 1995) LVG densities in eight sources have been cal- 
culated. Densities towards CS peaks within the 50" beam 
(0.56 pc at 2.3 kpc, the average distance of our sam- 
ple source) vary from source to source in the range: (3 - 
40)x 10 5 ctrT 3 . In four sources, the density falls by about 
an order of magnitude at ~0.8-2 pc from the CS peaks 
which could be connected either with the existence of lower 
density clumps or density gradients within clumps. The 
masses calculated from LVG densities are higher than virial 
masses and masses derived from continuum data implying 
small-scale dumpiness of the cores. The CS abundances 
have been calculated towards IRAS positions within 50" 
Gaussian beam for eight sources: X(CS)=(0.3-2.7)x 10~ 9 . 
The N 2 H + abundances for the same positions lie in the 
range: X(N 2 H + )=(0.3 -4.4)x 10~ 10 . 

3. For most of the objects the CS and continuum peaks are 
close to the IRAS point source positions. The CS(5-4) inte- 
grated intensities correlate with continuum fluxes per beam 
almost linearly in all cases while significant correlations be- 
tween the CS(5-4) and the N 2 H + (l-0) integrated intensi- 
ties are found only in five sources. The highest CS to N 2 H + 
integrated intensity ratios (~ 10) are found in the vicinity of 
the most luminous IRAS sources. The lowest ratios (S; 0.1) 
are observed towards positions close to the centers of the 
N 2 H + clumps without CS counterparts and IRAS sources. 

4. The study of spatial variations of molecular integrated in- 
tensity ratios to continuum fluxes reveals that for most 
of the sources the 7(N 2 H + )/Fi 2 ratios drop towards CS 
and dust emission peaks near IRAS sources, likely im- 
plying underabundance of N 2 H + towards these peaks. The 
/(CS)//(N 2 H + ) ratios are clearly enhanced towards the 
three most luminous sources, whereas the /(N 2 H + )/Fi. 2 
ratios drop towards these sources. Possible explanations 
of these results could involve dissociative recombination 
of N 2 H + or accelerated collapse processes. For CS(5^1) 
the 7(CS)/Fi. 2 ratios show no clear trends with distance 
from the CS peaks while for CS(2-1) such ratios drop to- 
wards these peaks implying line saturation effects in several 
sources. 

5. The analysis of normalized velocity differences between 
CS and N 2 H + lines has not revealed indications of system- 
atic motions in the sources towards CS peaks. A possible 
indication for infall motions has been found towards one 
position in G285.26. The CS(5^1) line widths in several 
sources have a tendency to decrease with distance from CS 
peaks which could be connected both with optical depth ef- 
fects and with enhanced dynamical activity towards centers 
of the cores. 
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Appendix A: The CS and N 2 H + LVG modeling 
results 

Using the LVG model we have examined the CS and N2H + 
excitation for a wide density range at distinct kinetic temper- 
atures. When calculating the excitation of the N 2 H + unsplit- 
ted rotation lines, we used the N2H + -He collisional rates from 
Daniel et al. (2005) (http://www.obspm.fr/basecoli. In Fig. IA.ll 
the CS(2-1), CS(5-4) and N 2 H + (l-0) brightness temperatures 
divided by N/V (the ratio of gas column density to velocity 
at the boundary) are shown versus density at two values of ki- 
netic temperature: 20 K and 40 K. For convenience of presen- 
tation they are multiplied by the factor of 10 14 and are denoted 
as P. They can be compared with observed ratios of molecu- 
lar line intensities to continuum fluxes (the latter are propor- 
tional to gas column density). The dependences of line optical 
depths (t) on density are also shown. We set d = 0.5 pc and 
AV = 2.5 km/s as typical clump size and line width (Paper I) 
and X(N 2 H + ) = 5 x 10" 11 and X(CS) = 10 -9 (which are close to 
the values for G285.26 from Table|6). These parameters enter 
the model as N/V = n(H 2 ) (2dX)/AV. 

The qualitative behavior of P versus m(H 2 ) is similar for all 
the lines, yet, for CS(5— 4) P reaches the maximum at higher 
densities than for CS(2-1) and N 2 H + (l-0). The position of this 
maximum depends on X(CS) (the lower X(CS) the higher is the 
density at which P reaches maximum). After the maximum, 
P decreases with density, approaching linear dependences (in 
logarithmic scale) due to line saturation effect. As the CS(5^1) 
lines in our observations are not saturated, we consider only 
densities 10 6 cirT 3 . Within the density range 10 5 — 10 6 cm -3 
P falls by a factor of <, 2 for N 2 H + (l-0) and by a factor ~ 2.5-4 




10 1 1 ^ 1 

10 4 10 5 10 6 10 7 10 8 

Q 

n(H 2 ) (cm ') 

Fig. A.l. The results of LVG modeling: CS(2-1), CS(5-4) and 
N 2 H + (l-0) brightness temperatures multiplied by the factor 
10 14 (A7 V)~' (P) versus density at two values of kinetic temper- 
ature. The dependences of line optical depths (t) versus density 
are also shown 



for CS(2-1) (Fig. IA.ll lower and upper panels). For the CS(2- 
1) lines such drop could be connected with saturation. 

We have no data on spatial distributions of kinetic temper- 
atures in our sources, yet, one-point CO(1-0) (Zinchenko et 
al. 1995) and IRAS data probably do not imply temperatures 
significantly higher than 40 K towards CS and dust peaks for 
the regions with angular sizes ~ 1'. If kinetic temperature rises 
from 20 K to 40 K towards density peaks, P decreases slower 
with the resulting decreasing factor of 1.3. Note that dust tem- 
perature gradients from 20 K to 40 K can reduce the observed 
7(N 2 H + )/F ratios by 2.3 times. Thus, one could expect that 
within the density range 10 5 — 10 6 cm -3 the 7(N 2 H + )/F ratios 
could drop by a factor not more than 3 as Tkin and 7^ rise from 
20 K to 40 K. 
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Appendix B: Infrared pumping and excitation of 

N 2 H+ 

Carroll & Goldsmith (1981) derived the following criterion for 
infrared pumping to have an influence on molecular rotation 
spectra: 

where / is a dilution factor of the pumping source which radi- 
ates at the temperature T s , A rot and A v ibr are the Einstein spon- 
taneous rates for rotational and vibrational transitions, respec- 
tively, v is the frequency of the vibrational transition, h and 
k is the Plank and Boltzmann constants, respectively. Taking 
f&R 2 /4r 2 , where R is the radius of the dust source, r is the 
distance from the source, and T S (K) » 51 (2 x 10 17 cm) 1/3 /T 1/3 
(Scoville & Kwan, 1976) one obtains: 

* [A~^ ( hv 1 R 1 ' 3 \ 
r ~ 2 V A rot 6XP 1 k 102 (2X10 17 ) 1 / 3 ]- 

The right-hand of this expression peaks at some R MAX , corre- 
sponding to: 

r MA x(cm) » 2.87 x 10 24 exp(-3) (y )" 3 , 

so that the pumping could be effective at distances not 
higher than this value. For the N2FF molecule with 
,u = 3.4D and v ro ,(l-0) = 93.174 GHz (Pickett et al. 1998) 
A rot » 3.6 x 10~ 5 s . For the first vibrational stretching mode 
with V!=3234 cm" 1 and A vibr = 857 s" 1 (Botshwina, 1984) 
r MAX ~6.9 x 10 15 cm «2x 10~ 3 pc. For the second bending 
mode with v^ = 698.6 crrT 1 (Owrutsky et al. 1986), the lifetime 
of the excited state is r = 123.1 ms (Heningeret al. 2003) corre- 
sponding to A v ibr = 8.1 s and leading to tmax ~ 6.6 x 10 16 cm 
« 2xl0~ 2 pc. For the third stretching mode with V3 = 2254 cirT 1 
(Botshwina 1984), r = 530.1 ms (Heninger et al. 2003) corre- 
sponding to A v ibr = 1 -9 s _1 and leading to r M AX ~ 9.5 X 10 14 cm 
~ 3 x 10~ 4 pc. 

A possible radiative pumping via rotational transitions due 
to the background radiation field, including IR source, is effec- 
tive only when the source optical depth at these frequencies is 
higher than unity. From our continuum data it is easy to esti- 
mate the peak optical depth at 1.2 mm averaged over the tele- 
scope beam in the optically thin case, 7 = F Q. 1 B v (JY l , where 
F is the peak flux per beam, £2 is the beam solid angle, B V (T) 
is the Plank function at frequency v and temperature T. Thus, 
for the strongest continuum source of our sample, G291.27, 
with peak flux F = 18.8 Jy beam 1 and T = 25 K, t * 0.05 av- 
eraged over a region with a linear size d~Q3 pc. Taking a 
power-law index of the column density-radius dependence of 
-1 which corresponds to the density-radius power-law index 
of -2, one obtains that in G291.27 the size of the region, where 
mean optical depth at 1.2 mm exceeds unity, is d~ 0.015 pc. 
Without averaging, the size of the optically thick region should 
be lower. Our calculations of CS and N2H + excitation assum- 
ing that infrared source radiates as blackbody with temperature 
Tir = 100 K (close on average to the hot dust temperature 



component of the sample sources, see Section l4~Tt . show that 
the dilution factor for the pumping source should be / <; 10~ 3 
in order to produce noticeable effect on the level populations. 
Therefore, the radius of affected region in G291.27 at 1.2 mm 
is r = R/(2yff) = d/(4yff) < 0.1 pc. The frequency of the 
N2H + (2-l) transition which mainly affects the J -I level pop- 
ulation is about 1 .5 times lower, therefore, the optical depth and 
the radius of the affected region should be 1 .5-2 times smaller 
(depending on the power-law index of the dust optical depth- 
frequency dependence). For the other sources this radius should 
be even smaller due to lower optical depths and lower distances 
in some cases. 



